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Based on a recent work by Schuller et al., a geometric representation of all skewonless, non- 
birefringent linear media is obtained. The derived constitutive law is based on a "core", encoding 
the optical metric up to a constant. All further corrections are provided by two (anti-)selfdual 
bivectors, and an "axion". The bivectors are found to vanish if the optical metric has signature 
(3,1) - that is, if the Fresnel equation is hyperbolic. We propose applications of this result in the 
context of transformation optics and premetric electrodynamics. 



I. INTRODUCTION 

In recent years there has been an increasing effort 
towards understanding and developing non-birefringent, 
exotic media. The motivation behind this trend appears 
to be twofold; on the one hand, many modern meta- 
materials - particularly in the context of transformation 
optics P, ~ are (bi-)anisotropic and yet strictly non- 
birefringent. On the other hand, the derivation of non- 
birefringent constitutive relations is a cutting-edge topic 
in the premetric electrodynamics of spacetime [3j-|5[ . 

Correspondingly, this paper finds its first application 
in the design of artificial media whose optical response 
is independent of the field polarization. Most notably, 
we derive a simple, relativistic constitutive law which 
parametrizes all non-birefringent, skewonless materials 
in terms of a symmetric matrix, two (anti-)selfdual 
bivectors and an axiorQ. Conveniently, the pivotal 
symmetric matrix is proportional to the optical metric, 
while the (anti-)selfdual bivectors vanish if the Fresnel 
equation is hyperbolic - that is, Lorentzian (3,1). 
By virtue of this result we also investigate on which 
transformation media can be attained whilst avoiding 
birefringence. Remarkably, it emerges that the key 
choice is to select vacuum as the space on to which 
to perform the ray-guiding coordinate changeQ. Other 
options are recognized to be viable. However, they 
require introducing an axion and/or a metric signature 
different from (3, 1). 
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A "skewon" arises, for example, when the permittivity tensor 
or the permeability tensor puj are not symmetric. An "axion" , in- 
stead, contributes to the medium response with a non-reciprocal, 
isotropic term. See Section [H] for further details. 
The vacuum ansatz can correspond to a curved spacetime - pro- 
vided the appropriate refinements are considered [f|. 



A second application of this paper targets the empty 
space structure of electrodynamics, rather than exotic 
materials. The seminal works on this topic - due to 
Hehl, Obukhov and Rubilar 0, 0, ls| - demonstrated that 
the light-cone structure of spacetime can be re-derived 
using some minimal, premetric, experimentally justified 
assumptions. In particular, one obtains the conventional 
vacuum response merely by setting the skewon part to 
zero, while enforcing a specific closure condition. A 
similar, yet separate,scheme is pursued by Lammerzahl, 
Hehl [4j and Itin (5|. They achieve the customary 
light-cone by ruling out birefringence and requiring that 
running-wave solutions must exist in all directions. So 
far, no link has been made between this latter approach 
and the former skewonless, closure-abiding one. Our 
work presents a direct connection, at least for media 
with no skewon. 

This paper is organized as follows: firstly we introduce 
bi-anisotropic local linear materials, the Fresnel equa- 
tion and the appropriate relativistic toolbox (SectionllT]). 
Subsequently, a classification of all non-birefringent me- 
dia is obtained and the corresponding compact represen- 
tation is promptly developed (Sections MI AHhTCI and 
the worked example of Section Till D[) . The results are 
then applied both in the context of premetric electrody- 
namics (Section IIV A[) and in the context of transforma- 
tion optics (Section llVBp . We draw our conclusions in 
Section [V] The Appendices lATICl summarize some tech- 
nical aspects, mostly concerning the work of Ref. @. 



II. FRESNEL EQUATION AND 
NON-BIREFRINGENCE OF LINEAR MEDIA 

Since this paper is partly targeted to applications in 
engineering, we assume the existence of a "background" 
metric g, thus confining all premetric refinements to Sec- 
tion IIV Al We consider linear, non-dissipative media, 
for which we introduce a relativistic (frame independent) 
representation. Suitable starting points are the canonical 
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Boys-Post relations 
D =T E B 



H = /3 ■ E + Jl~ x ■ B , (1) 



which define the permittivity e, the inverse permeabil- 
ity /I -1 and the magneto-electric couplings a and /3 [Icj . 
More compactly, one can also encode the medium re- 
sponse asO 



(2) 



The 6x6 block matrix, which is identified here, is often 
labeled x IJ ~ with / and J varying between 1 and 6. Fur- 
thermore, it is common practice to introduce the identi- 
fications (Fj) = (-E; B) and (W 1 ) = (D; H). This 
allows one to re-express the map © as0 



D 




—e a 




—E 


H 








B 



w 1 



X 



I J 



F 



J ■ 



(3) 



A relativistic notation is now within reach and is achieved 
by implementing spacetime indices, ranging from to 
3 and commonly denoted by Greek letterqj. With this 
notation in mind, Eq. ([3]) is immediately translated to 



(4) 



where the quantities F a p and W a @ are each defined 
through an appropriate tableau [Tl| : 



F 



a/3 



— Ei — E2 — E% 

Ei B3 —B2 

E2 — -B3 Bi 

E3 B2 —B± . 

■ D x D 2 D 3 

-Di H 3 -H 3 

-D 2 -H 3 Hi 

-D a H 2 -Hi . 



(5) 



(6) 



By inspection - or by observing that F a p and W a @ are 
antisymmetric - one can conclude that each tableau con- 
tains only 6 independent entries. These can in turn be 
collected according to the rule 

{[01], [02], [03], [23], [31], [12]} -j. {/ = 1, 2, . . . , 6} , (7) 

so as to retrieve the 6-dimensional vectors Fj and W 1 . 
Using an analogous argument, since 



one can always switch between the spacetime form ^ va P 
and the 6x6 form \ = {y/ J }E1 As one might expect, the 
relativistic representation is particularly convenient when 
considering a change of frame. More specifically, given an 
arbitrary non-singular transformation matrix L a a , one 
has that (cf. Sec. MM 



X 



fj,' ' U' ' OL '/?' 



t ft' t u' j a' j 0' fivafi 
L fi L v L /3 X 



(9) 



and that the tensors F a p and W a ^ behave in a similar 
manner. For many applications, it is also useful to de- 
compose x by means of the matrix identity 



= x+x T . x~x T 

X = n 1 o 



(10) 



whereby one can isolate a symmetric contribution x Sjmm 
and an antisymmetric contribution ^x Q- A further 
split, reading 



Xs 



(D f+ (3) = 



(11) 



can be achieved by requiring that W^l/"" 3 ^] vanishes, 
while ( 3 ) x ltM/al3] does nofl By virtue of HTDJ) and (HJ), one 
can finally separate the medium response into a "princi- 
pal" part, a "skewon" part and an "axion" part: 

ylivafS _ {l)^ftuaf3 _|_ {2)^fiuaf3 ^_ {3)^ftvaf3 ^2) 

where ( 3 ) x ^ a £ = a& va() , and ^ va ^ is the Levi-Civita 
tensor defined via 



fivotfS 



£(3123 



■det(f)] 



(13) 
(14) 



A non-zero ^x component is invariably excluded, since 
it implies that any of the following common symmetries 
is broken: 



e = e T 



a=-/3 



(15) 



By contrast, a finite ( 3 ) y is still quite rare, but has been 
observed in nature [lU, [l3[ - thus overturning a popu- 
lar dogma ("Post's constraint", see dUdi]]). The axion 
field also finds a practical application in the perfect elec- 
tromagnetic conductor (PEMC) proposed in [l5[. Given 
these considerations, the present paper assumes that all 
skewon contributions vanish. 



^luotf) _ _yVfia(3 _ _^ltvf3a 



(8) 



3 The minus sign appearing in front of E is clearly superfluous - 
but only for the purpose of this equation. In many other occa- 
sions it is required so as to comply with Lenz's law 

4 Throughout this paper we make use of Einstein's summation 
convention. 

5 Given a coordinate patch time is parametrized by x°, while 
space is spanned by the components x 1 , i = 1, 2, 3. 



6 Additional information about this topic can be found in [Tl| . 
In all that follows, we use the relativistic and 6x6 notations 
interchangeably. 

7 The reason for labeling the antisymmetric contribution ( 2 ) will 
become clear soon. Ultimately, we are aiming to match the no- 
tation of 

8 The square brackets used here denote index alternation [Till . One 
should also notice that the permutation 

is identically zero. 
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The easiest way to specify the properties of a mate- 
rial is via the 36 entries of x- However, when a symbolic 
calculation must be carried out, it is often more conve- 
nient to reduce the constitutive relation to an abstract, 
more compact form. In engineering, this is mostly at- 
tained by simplifying e, ~p, a and (3 separately, according 
to their individual symmetry properties [16| . This strat- 
egy, however, is manifestly not relativistic, as it involves 
manipulating the four 3x3 sub-matrices of ^ indepen- 
dently. Therefore, at least for the purpose of this work, 
one must choose a different approach, where 
composed in terms of simpler spacetime quantities. For 
instance, vacuum is characterized by 



(16) 



where g a P is the inverse of the background metric g a p. 
Relaxing some restrictions leads to a first class of non- 
birefringent materials: the Q-media [l7| . These are given 
by 



X 



Q — 



(17) 



with sq = ±1, and provide three supplementary sources 
of flexibility. Firstly, the tensor Q a " is not a metric, 
since although it is invertible, it is not necessarily sym- 
metric. In the context of this paper, this property has 
little relevance, since a vanishing skewon implies that ei- 
ther Q = —Q T or Q = Q T . Insomuch as thejformer case 
leads to an unspecified wave (co-) vector |l8(, only the 

latter case can be accepted. Thus, Q can henceforth be 
considered to be a "constitutive" metric, converging to g 
for the vacuum (TIB")) . Accordingly, media of type (TTf]) will 
be referred to as "metric media" . Secondly, the consti- 
tutive signature need not be (3, 1), that is, "Lorentzian" . 
This attribute will prove essential in deriving a compact 
representation for all non-birefringent materials (Section 
Mil) . Lastly, the impedance of the medium in Eq. (|17[) 
can be calculated to be 



|det(Q^)| 
-det(g Q/3 ) 



(18) 



and is not necessarily equal to (/loAo) 2 - Moreover, the 
siffl i sq allows one to obtain a negative index of refraction 

A further class of materials, which naturally extends 
Q-media, has been proposed by Wallen and Lindell 
[l3,[23|. They modify (HZJ) by including two "bivectors" 
(antisymmetric tensors) A and A: 

xZZ 13 = s q (Q^ a Q uf) - Q^Q ua ) + A^A a ? , (19) 

where, in this article, one must enforce A — saA with 
sa = ±1) since ^x i s assumed to vanish. 

Formulating the constitutive relation is the first step in 
characterizing the behaviour of light. The second mile- 
stone consists in solving Maxwell's equations - in general 



or, as in the present work, in two specific circumstances. 
When studying the structure of exotic spacetimes, the 
propagation of sharp electromagnetic fronts is considered 
[3j. By contrast, when modeling laboratory materials, 
some approximate plane-wave analysis is implemented 
(the light's front-velocity in a table-top medium is always 
exactly c [2ILI22I]). Both of these scenarios correspond to 
the so-called "geometrical optics" limit, whence a disper- 
sion ("Fresnel") relation, linking the angular frequency 
lj to the spatial frequency k, is derived. Using frame- 
independent notation, the 4-wave co-vector K = (— uj, k) 
is found to obey 



a ^ 6 ^^ K X 5prTX K„.K v K n K n 



0, (20) 



that is, a quartic Fresne l poly nomial, cubic in the 
medium parameters [!, H, l23T - [25l | . By using Eq. (TT2"|) . 
one can also prove that the effect of the axion on this 
formula is zero. Consequently, it impossible to specify 
^X by analyzing light propagation in bulk. 

A material is non-birefringent if the Fresnel polynomial 
becomes bi-quadratic: 



(G a P K a Kp)' z = 



(21) 



where G is an "optical" metric. For an explicit deriva- 
tion of the non-birefringence conditions we refer to the 
literature [4, 5] (cf. Section UTTDj). The Fresnel equation 
for a Q-medium is trivial (G a/3 ocQ al3 ). However, that of 
its extension (fl9)) is not [10(: 



Q — saA ■ Q~ ■ A) K a K p 



a/3 



= . (22) 



Birefringence, as induced by the bivector A, can be elim- 
inated by requiring that 



A 



A = aQ , 



(23) 



where the proportionality constant a may vanish, but 
a 7^ sa- The generalized Q-media that satisfy (l23|) cover 
only a range of non-birefringent materials. Nonetheless, 
they provide essential guidance for the generic case. 



III. GEOMETRIC REPRESENTATION OF ALL 
NON-BIREFRINGENT MEDIA 

A. From Schuller's classification to a compact form 

Discovering specific examples of non-birefringent me- 
dia is relatively easy. However, finding a general repre- 
sentation is - at a first sight - difficult, given the com- 
plicated structure of the Fresnel equation ([2"U|) . Re- 
markably, a classification scheme recently proposed by 
Schuller et al. brings this task to the realm of possibil- 
ity (at least for a vanishing skewon field) . The key break- 
through comes from an in-depth analysis of the Segre 
types |26l . l27l ] emerging from the modified tensor 



per 



1 



(24) 
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As an outcome, 23 equivalence classes are obtained, 
each equipped with a normal form (or "typical repre- 
sentative"), whereby a complete categorization is readily 
achieved - an arbitrary "symmetric" constitutive ten- 
sor can be linked uniquely to a typical representative 
by using a change of frame. Simultaneously, the non- 
birefringence conditions developed in 0, Q are invariant 
under any basis transformation, as dictated by the princi- 
ple of relativity. Therefore, an observer independent cov- 
erage of all non-birefringent media can be fully attained 
from the 23 specific matrices of Ref. _9j - complying, in 
this sense, with the main aim of this section. 

A few technical comments are in order. Schuller et al. 
opt to classify the lower- indexed "area metric" Gp Upai 
rather than the upper-indexed x fJ,1,pcr which is favored 
here. Regardless of this choice, the normal forms that 
they derive do apply directly to \. In addition, the 
representatives in Ref. [9| are subject to the condition 
detx = 1- We choose to temporarily suspend this con- 
straint, since the procedures below never preclude us 
from achieving a unit determinant. As a final remark, 
one should notice that Schuller et al. employ an alterna- 
tive 6x6 index ordering - compare equation (2) in Ref. 
Q with equation (J7J of this manuscript. 

With the help of computational software^, the re- 
quired non-birefringent classification is derived straight- 
forwardly. The immediate outcome is a set of matrices, 
which provide a valid solution, but not a transparent an- 
swer (see Appendices lAl and iBj) . In spite of this difficulty, 
one necessary - but not sufficient - condition does emerge 
clearly. The structure 



= s Q (Q w Q va - Q pa Q up ) 

+ s A A pv A pa + s A A pu A p,J + ae^f" 7 , 



(25) 



encompasses efficiently all skewonless non-birefringent 
materials, although some polarization-dependent solu- 
tions are yet to be ruled out. A complete refinement will 
soon be obtained. Nonetheless, some essential remarks 
should not be avoided. Equation (f2"5"j) is designed to con- 
veniently encode the optical metric, solely by means of 
Q. Moreover, our compact form naturally generalizes Eq. 
(|19[) . at the small price of including an innocuous axion 
offset. It is also important to realize that the representa- 
tion (|2"5|) is invariant under some "bivector space" trans- 
formations comprising A and A. If sa s I = 1) trigono- 
metric rotations of the form 



cos< 
sin< 



— sin i 

cos <A 



(26) 



map the medium tensor onto itself. If s^s^ = —I, proper 
rotations have to be replaced by hyperbolic rotations 



A'\ _ ( cosh ip — sinh ip 
A' I \ — sinh cp cosily 



(27) 



B. Necessary and sufficient conditions imply duality 

At this point, the crucial idea is to study the Fresnel 
equation of x^be", so as to find a true, minimal, non- 
birefringent representation. From Eq. (|20| : 



det[Q]((Q - s Q s A A ■ Q- 1 ■ A) pv KpK v y. 



x (Q-s Q s A A 



sas a 



r 1 • A) pa KpK a 

n 2\ 



(A • Q" 1 • Af v KpK v = . (28) 



Thus, the medium is non-birefringent if and only if the 
following conditions are met: 

A-^- 1 ■ A-^- 1 = ai t , (29) 
A ■ T}^ 1 ■ A ■ = a 2 t , (30) 
A ■ (J" 1 ■ A ■ "Q- 1 + A ■ (J" 1 • A ■ T^- 1 = 2a 3 l . (31) 

Here, the scalars a\, a 2 and 03 must always comply with: 

(1 - s Q s A ai){l - s Q 8 A a 2 ) - s A s A a 3 ^ . (32) 

A solution to the above constraints is easily found by 
exploiting the bivector identity 



detQ| 

where eZ ' g is a Levi-Civita-type tensor, with dual e 



e?i2 3 = |det(Q)|- 



( £ Q)^ = sgn[det(Q)](e Q 1 ) /J 



(33) 

3 : 

(34) 
(35) 
(36) 



Applying the relation ^ in Eqs. (gU) and (SOI), yields 
our main results; namely, 



A pu _ fX_ el J.uo .',77-! A — 



{Q^-A-Q~ V U -=s x (*AY v 



(37) 



A"" 



S Y pvufl 

2 e Q 



(Q- 1 ■A-Q- 1 ) ar .= s Y CAr» , (38) 



that is, A and A must be selfdual (sx, sy = 1) or anti- 
selfdual (sx,sy = — 1), with: 



A = **A , 



A = **A 



(39) 



The additional condition (|3Tj) is solved automatically, 
provided sx = sy- 

Rotating the constraints (I2^|) - (|5T|) by means of Eq. 
(|26p . determines that 03 = can be achieved whenever 
sas a — 1. By contrast, if sas a = —I, the hyperbolic 
rotations (|27[) allow one to attain a 3 = only when 



The results of this paper were calculated using Mathematical 



ai + a 2 
03 



i [-2, +2] 



(40) 
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C. A summary with all explicit solutions to duality 

In summary it is found that the optical metric of any 
non-birefringent medium obeys cx Q^ v . All non- 
birefringent, skewonless materials are characterized by 
a Q-medium "core" , which encodes the optical metric. 
Two (anti-)selfdual bivectors and an axion term can be 
included as corrections, but they do not affect the disper- 
sion relation. Any further discussion needs to distinguish 

between the three possible signatures of Q. 

1. Signature (3,1) 

When the signature of Q is Lorentzian (3,1), the fol- 
lowing identity is true for any bivector xj) [3j: 

V> = . (41) 

Correspondingly, equation (j39|) cannot be satisfied, un- 
less A = and A = 0. Then, any non-birefringent mate- 
rial, with hyperbolic Fresnel equation, is simply a metric 
medium - complemented by an eventual axion. Some 
observations in this sense can also be found in the group 
theory analysis of Ref . [9J . 

2. Signatures (4,0) and (2,2) 

The remaining two signatures support finite (anti-) 
sclfdual bivector corrections. These can be written in 
terms of the vectors t, u and v, such that: 

A v> = 2rK^ + sxego*®- 1 ■ x) p (^ 1 • u) a , (42) 
ft" = 2 >^1 + sxe^^- 1 ■ x) p (^- 1 • v) a . (43) 

In this parametrization, r singles out a particular direc- 
tion, according to which all other vectors can be decom- 
posed. In this sense, r plays a role similar to that of a 
time vector, and can be selected independently of A and 
A - in almost all cases. Complications arise only if the 
signature is (2,2) and r is light-like, a scenario that we 
shall resolve later. Thus, one can take 

t-Z^-t^O, (44) 

while no specific signature is assumed yet. The compo- 
nents of u and v parallel to r never contribute to Eqs. 
(l4"2l and P3")l . Therefore, given that we excluded a light- 
like r, one can choose r • Q _1 ■ u = and r • Q^ 1 ■ v = 0, 
with no loss of generality. With these simplifications, it 
is possible to derive the following relations: 

ai=-(T- Z Q- 1 -T)(u- : Q- 1 -u) , (45) 
a 2 = -(t -t}- 1 ■ t){v -l^- 1 ■ v) , (46) 
a z = -{T- Z Q- l -T){u- Z Q- 1 -v) . (47) 

If Q has signature (4,0), a true norm is established. 
Thus, the occurrence of a\ — or a 2 — is forbidden, 
unless t, u or v vanish. By means of trigonometric or 



hyperbolic rotations, the orthogonality u -Q^ 1 v — can 
be obtained without altering x NBR . Equivalently, = 
can always be assumed. Hence, for signature (4,0), the 
three vectors r, u and v can be chosen to be mutually 
Q-orthogonal, with no effect on the medium response. 

If Q has signature (2,2), the situation is more compli- 
cated. The vectors it and v can be light-like, thereby 
allowing a\ — and a 2 = 0. Moreover, if sas^ = — 1, 
one can find cases where 03 = is no longer achievable. 

Finally, we should comment on the possibility of a 
light-like r. This choice is more subtle since there ex- 
ist pairs of vectors (t, u) obeying all of the following 
conditions: 

>u^0, t-^- 1 -u = 0, (48) 

t- Z Q- 1 -t = Q, u Z Q- l u = Q. (49) 

Pairs of this kind span the kernel of the map (|4"2"j) , which 
links two vectors to one bivector. Whilst in specific sit- 
uations one could allow for a light-like r, this choice is 
inconvenient in general. 

D. A worked example: revision and further details 

A worked example is proposed, which provides hands- 
on revision for Sections |III Al InTCl By implementing the 
prescriptions of Refs. [U, a specific "Schuller" normal 
form is explicitly reduced to the non-birefringent limit 
(1231) . Two bivectors, playing the role of A a ° and A a @ , 
are readily determined. Their structure is shown to be 
encompassed by Eqs. (|42|) and (|43|) . thus illustrating a 
specific instance of (anti-)selfduality. 

As explained in Section IlIIAl choosing a typical repre- 
sentative singles out a Segre type. In turn, this identifies 
a set of eigenvalues for the modified medium tensor of 
Eq. (|2~4")l : in this example, we select the characteristics 

o-i+in, 0-1 -in, Ai, A 2 , A 3 , A 4 , (50) 

where a%, n, Ai, . . . , A4 are real scalars. The associated 
6x6 matrix normal form Xnf reads 

" -n (7i " 

A 3 A 4 

- = Ai A 2 

Xnf (7i n ' { ' 

A 4 A 3 
.0 A 2 Ai . 

as reported in Ref. [9] (cf. "metaclass VI" therein, with 
the index ordering ([7|l). At this point, the optical prop- 
erties of the material are fully specified and the Fresnel 
equation (|20[) can be invoked to obtain the polynomial 

M q 4 + Ahq 3 + M 2 q 2 + M 3 q + M 4 = , (52) 

where q is a component picked arbitrarily from the wave 
co- vector K^. Crucially, the coefficients Mo, . . . , M4 are 
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Quartic Equation 



M =0 




M =0 
& 

M=M 2 /(4M.) 



Y 



M 4 =(4M M 2 -M 1 2 ) 2 /(64M 3 ) 
& 

M 3 =M 1 (4M M 2 -M, 2 )/(8M 2 ) 



N 



N 



Non birefringent media | All other media 



FIG. 1. The procedure to make Eq. (|52p a bi-quadratic. One 
must check that A/4 7^ at every node, so as to avoid q = 0. 
The birefringence elimination scheme proposed in Refs. [3, [^| 
does not include the Mo = branch depicted on the left. 



still a function of (— w, fc); however, they are independent 
of the particular entry K u such that q — K u . For the 
material under consideration, if q coincides with K\ : 



M = - A1A3T1 , 
Mi = , 

M 2 = + Kin [-A 2 + (A 2 - A 4 ) 2 - A 2 ] 

,2 , T 2" 



(53) 
(54) 



,2 , T 2 



U 

f] • 



(55) 



-A 3 2 A 3 [-A 2 + (A 2 - ( t 1 ) 
-Jf|Ai [-A 2 + (A 4 - <ti) 
M 3 = - 2K Q K 2 K 3 x 

{+A 2 (A 4 -a 1 )-A 2 (A 4 -a 1 ) 

- a\ [A 2 - A 4 (A 4 - at)] 

- A 4 r 2 + A 2 (A 2 - A 2 + a\ + r 2 ) } , (56) 
M i = -K$\ 1 \ 3 T 1 + X 2 A1A3T1 +AtA 1 A 3 r 1 

+ A 2 2 A 3 2 n [A 2 + A 2 - (A 2 - A 4 ) 2 ] 
-A 2 A 2 2 A 3 [ r 2_A 2 + (A 2 -a 1 ) 2 ] 
-X 2 A 3 2 A 1 [r 2 -A 2 + (A 4 -a 1 ) 2 ] , (57) 

as one can verify by using a numerical software (see Foot- 
note EJ). Two topics ought to be addressed: firstly, the 
identifications ([53]) -([5*7 ]) depend on the choice of q, in 
contrast with the relativistic standard retained so far; 
secondly, birefringence must still be eliminated through 
the reduction of Eq. ([55) to the form (l2"Ti 

Remarkably, once a bi-quadratic (G Q/9 K a Kp) 2 = is 
derived, the result is independent of how q is defined. 
This fact can be verified directly, by means of a computer, 
or geometrically from the uniqueness of the "light-cone" 
Knowing that (when birefringence is ruled out) a 
frame independent outcome is ensured, one can safely 
apply the scheme of Figure [JJ and retrieve Eq. (PHI) . 

The analysis underpinning this method is due to 
Lammerzahl, Hehl Q and Rin The present exam- 
ple demonstrates that, in extension to the existing the- 



ory, non-birefringent solutions can be obtained even for 
Mo = (Fig. [JJ left-hand branch). By contrast, M 4 must 
always be nonzero, otherwise the dispersion relation ([52]) 
leads to the unphysical root q — 0. Carrying out the pro- 
cedure described in Figure [JJ for the medium ([51]) yields 



M 


= 0, 


Ai 


= 0, 


A 2 


= ±A 3 - 


I-A4, 


(58) 


M 


-0, 


A3 


= 0, 


A 4 


= ±Ai - 


f-Aa, 


(59) 


M 


= 0, 


Tl 


= 0, 


ci 


= ±A 3 - 


h A 4 , 


(60) 


M 


-0, 


Tl 


= 0, 


0"! 


= ±A X - 


hA 2 , 


(61) 



whilst taking Mq 7^ provides no results compatible with 

(4M M 2 -M 2 ) 2 
+ 1 °J a , 3 ^ ■ (62) 



64M 3 

It is worth noticing that the eigenvalues in Eq. (|50[) can be 
re-ordered, for example by exchanging A 3 with Ai and A 4 
with A 2 . Accordingly, one can observe that Eqs. (f5"9"]) and 
(IBTT) are equivalent to Eqs. (|58p and (l60l) . respectively. 
With no loss of generality, it is thus possible to focus on 
the two cases 



(i) 
(ii) 



A 3 =0, 
ri =0, 



A 4 — — £Ai + A 2 
<J\ = +CAi + A 2 



(63) 
(64) 



with C = ±l. The corresponding Fresnel equations ([20]) 
are non-birefringent, as required, and read 

-Ai[(A 2 -a x - CA1) 2 + r 2 ](AiA 2 + CA A 3 ) 2 =0, (65) 
-Ai[(A 2 - A 4 + CA1) 2 - A 2 ](AiA 2 + CA A 3 ) 2 =0. (66) 

Moreover, the pre-factors multiplying (K\ K2 + CK0K3) 2 
must be non-zero, otherwise the dispersion relation van- 
ishes identically. One thus enforces the conditions 

(i) A 17 ^0, (A 2 -ai-CA!) 2 +T 2 ^0, (67) 

(ii) A 17 ^0, (A 2 -A 4 + CAi) 2 -A 2 ^0, (68) 

so as to identify the optical metrics for cases (i) and (ii) : 



G (i) — G ( 



/o 1 
c 
c 

V 1 



(69) 



As an aside, one should observe that G (i) and G (ii) have 
signature (2, 2). Hence, these metrics are not suitable for 
governing the causal structure of spacetime [9, [28| . All 
realizations of Eq. ([6"9"]) are restricted to laboratory mate- 
rials, in the steady-state limit (Section [IT]). In agreement 
with the findings of Section IIII AI the non-birefringent 
media (f6"3")) and (f6"4"]) have a constitutive relation of the 
form ([23]). More specifically, one can show that 



A(i) 



A(ii) 



sgn(Ai) (Q^Q^f 



-sgn(Ai) {Q^QIZ-QZQl 



(70) 



+ s^A^AZ + hiAmAZ + x ^ Upa > ( 71 ) 
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where the two constitutive metrics <5m and Q"^ coincide can demonstrate that 

*L = VK\ , !L = VK\ s (u) , (72) 



and are non-singular by virtue of Eqs. (joTl) and (1551) . The 
bivector corrections and A {ii) are given by 



a w 

-O(0 







C (i) 

-C (i) 










-a (u) 





-c (ii , 








b m 






d w 

-d m 

o (il) 

b (m 





(73) 



(74) 







c (ll) 
d (u) 

-d(ii) 



(75) 
(76) 

where the real scalars {a (i) , . . . , d (i) } and {a (ii) , . . . , cZ (ii) } 
must be expressed in terms of the eigenvalues (|50p . Simi- 
larly, the signs {s (i) , S (i) , s (a) , S (ii) } are still to be specified. 
One can directly verify that, defining the variables 



= -C^i - A 2 +(7l , 
= +CAi - A 2 + A 4 , 
= sgn(|A 3 |-|p ( a)|) , 



(77) 
(78) 
(79) 



and selecting the signs s (i) , s (i) , s (ii) and s (u) according to 

s w =sgn(n), s (l) =- sgn(ri) , (80) 

s (iJ) = sgn(A 3 ) , s CJS) = C ( ,)Sgn(A 3 ) , (81) 

is compatible with the following choice of bivector entries: 



6(0 
c (i) 
d (i) 

C (tt ) 
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(ii) 



d. 



(") 



-d (i)) (82) 

+ c w , (83) 

+ sgn( Tl ) [(| P(1) | 2 + Inl 2 )^ + |n|]4/>/2 , (84) 

-8ga{p w )[(\p m \* + \Ti\*)i-\Ti\]i/y/2, (85) 

+ C (l i)Sgn(A 3 )sgn(p (li) )6 (il) , (86) 

-C(tt)Sgn(A 3 )sgn(/9 (11) )d (tt , , (87) 

+ [|A 3 |-e (li ,|p (li) |]VV2, (88) 

+ [Ki)l+C( U )|A 3 |]'/^- (89) 



As a final consistency check, it is possible to verify that 
the conditions (|3"?| and (|3"5)) are satisfied. In fact, one 



A UV 






T (i 


\ — 1 

))pW(i) ' 


u (i)J(X i 






+ c^7(^ 


T(i 


))p(Q(i) 1 










T (i 


i))p( ( 3(ii| 


•W(u)) CT 


K 


— ZT (ii)f(ii) 




T(i 


l))p(Q(U) 


■«(u))a 



where the triplets {T (i) , u (i) , and {t (u) , it (ii) , are 
specified (non-uniquely) as 

t w = (0,1,0,0) , (94) 

u w = (o w /15,0,0,6 w /17) , (95) 

t> w = (c (i) /15,0,0,d (i) /17) , (96) 

r (u) = (0,0,1,0) , (97) 

u (u) = (o (u) /17,0,0,V)/15) , (98) 

t> (tt , = (c (I1) /17,0 J 0,d (tt) /15) . (99) 

Clearly, the bivectors d73l - (|7p]) display the typical struc- 
ture (|4*2")l - (j4*3")) . Thus, the corrections A w , . . . , A Cil) are 
sclfdual for C = +1) and anti-selfdual for £ = —1. 



IV. UNIQUE LIGHT-CONE AND 
TRANSFORMATION OPTICS 

A. The light-cone and premetric electrodynamics 

This section is organized as follows: firstly, some re- 
finements are introduced, as required by the premetric 
electrodynamics of spacetime. Then, it is demonstrated 
that the results of Section IIIII remain valid, besides all 
amendments. Finally, the work of Ref. 0] is linked to 
that of Ref. 0. 

Following the premetric formalism of Post [11| . the 
mathematics used so far needs one adjustment. The per- 
mutation symbols e a ^^ s and e a p 7 s, with 



„012 3 



CQ12 3 — 1 



(100) 



must be employed. Correspondingly, it is necessary to 
go beyond the rules of tensor transformation, since 



e «'£V<5' = \det{L p ')\- 1 L a ^L f3 'L\L s ' s e a ^ s 



(101) 



where the "density" factor \det(L p p )\~ w is said to have 
"weight'' w — +l. Given that the premetric medium re- 
sponse x PM i s a derived quantity [3j , 



livuj) — \pV>vpv a/3 
Apm ■ 2 P a 



(102) 



the following "amended" transformation rule is obtained 
(compare with Eq. ([9])) 



Apm 



PM 



(103) 
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so that X™ is a tensor density with weight w = +1. 

At this point, we verify that the refinements (jlOlj) and 
(|103p do not affect the findings of Section Hill Crucially, 
one can prove (cf. Appendix [C]) that Schuller's classifi- 
cation still applies, and that the representation (f25j) is 
translated to: 



-V PM 



det(Q^)|M[(QWQ w - Q^ a Q vp ) 
saA^A"" + s A A pu A pa ] + ae^ pa , (104) 



where M is a true scalar. The "differences" between Eq. 
(I104p and the counterpart (12~51) drop out from the Fres- 
nel equation. Hence, the expression (|2"51) is always ob- 
tained, together with the conditions (157)) and ([55]) . Con- 
sequently, all the conclusions of Section UTT1 remain valid. 

To understand the importance of our work in the con- 
text of premetric electrodynamics, a few results of this 
program should be mentioned. Hehl and Obukhov Q 
consider some key experimental facts, and thus deduce a 
constitutive metric - a "light-cone" structure - emerging 
from the response x PM - By enforcing electric- magnetic 
reciprocity, they require the axion free contribution xHi 
to obey the closure condition 



e IKXpM e -LMXp M 



ki K kk J 



(l/Z 2 )5/, (105) 



where Z is a constant. As a second and final constraint, 
they also set the skewon part to zero, so as to obtain the 
desired Lorentzian metric - together with an arbitrary 
axion contribution. In Ref. Q Lammerzahl and Hehl ex- 
ploit different observational facts and develop a parallel, 
yet unrelated, derivation. Rather than simplifying the 
material response directly, they study the characteristics 
of the Fresnel equation, thus achieving an optical metric. 
They impose that: 

• The medium is non-birefringent, that is, there ex- 
ists a factorization of the form (|2~Tj) . Hence, an 
optical metric is uniquely defined. 

• Given a value of lo, the optical metric supports two 
real ( "running- wave" ) solutions - accordingly, the 
Fresnel polynomial is hyperbolic. 

Clearly, the two schemes presented above are conceptu- 
ally distant. Ref. [1] proceeds by exploiting symmetries, 

and deduces a constitutive metric Q. By contrast, Ref. 
|4j takes a strictly "eikonal" route, and derives an opti- 
cal metric G. Given these differences, it is remarkable 
that, in the absence of skewon, we can indeed establish 
a link: as seen in Section IIII1 the optical metric for all 
non-birefringent media is proportional to the constitutive 
metric, Q a/D cx G a/3 . Observing that vacuum must have 
signature (3, 1)E3, the result A = A = is found imme- 
diately. Thus, the distinction between constitutive and 



optical structures disappears: 

z- 1 



APM 



V- det G a t 



\G w G l " J -G >ia G vp \+ae llupa , (106) 



up to a choice of Z and a. These residual sources of 
flexibility do not distinguish the works of Ref. Q and 
Ref. [3] . Indeed, Hehl and Obukhov "manually" set Z = 
Weo) 1/2 and a = in (D.6.11) and (D.6.12). 



B. Non-birefringent transformation optics 

Transformation optics (TO) is a useful design tool, 
which offers direct control over the light rays [l], 0]. 
Its working principle is simple and can be explained 
in few sentences. The starting point is always a well- 
understood setup, called "base-geometry", supporting a 
pre-determined electromagnetic configuration. A coor- 
dinate change, or "deformation" , is then applied, which 
drags the fields to a useful layout - and which specifies 
the necessary medium parameters. Thereby, a powerful 
recipe is spelled out, which can be used to engineer ma- 
terials with tailored optical properties. Novel devices can 
thus be developed in the laboratory, whilst dissolving any 
reference to the enabling grid-distortion. 

Most applications of TO, such as cloaking [il l29l. l30l| . 
light- harvesting [31[ and lensing |32[ , make use of a sim- 
ple base-geometry: usually vacuum, or an arrangement 
of scalar e and fi. By virtue of this choice, birefringence 
is always ruled out from the very outset and never affects 
the design process. For many purposes, this is a true ben- 
efit, which should be exploited fully; whence the question: 
what is the most general substrate for constructing a TO 
scheme free of birefringence? 

In the case of zero skewon, a complete answer comes 
from the analysis of Section IIIII In particular, any ma- 
terial featured in the base-geometry must take the com- 
pact form ([251 , whilst obeying the duality rules (1571) and 
(1581 . A number of interesting scenarios are encompassed 
in this general solution - for example, when the signa- 
ture is (4,0) or (2,2). Yet it appears that, technolog- 
ically, the Lorentzian metric does enjoy a special sta- 
tus, since it guarantees propagating bulk modes. Hence, 
under most circumstances, the Fresnel equation is hy- 
perbolic, and the bivectors A and A vanish. Accord- 
ingly, non-birefringent TO is reduced to the simple class 
of (3,l)-metric media (1171) . complemented by an axion. 

Given this last consideration, it is easy to understand 
why avoiding birefringence requires the TO substrate to 
be vacuum, with only few exceptions. A comparison with 



In empty space, causality demands the Fresnel equation to be 



strictly hyperbolic [sl |28|| . For table-top materials, this restric- 
tion can be relaxed, since we assumed continuous wave operation 
- rather than information transport, as encoded in the front ve- 
locity (Section Hill. 
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Xo m Eq. (fT())) reveals that, at this point, only two quan- 
tities can still be adjusted: the impedance (TTHI) . with the 
sign sq, and an axion term. Even so, when probing an 
interface with empty space, selecting 



/det(Q^) 
det(5 Q/3 ) 



(107) 



or choosing ^x 0j generates finite reflections in trivial 
scenarios (see, for example, Ref. [33]). By observing that 
TO does not naturally account for back-scattered waves, 
one can conclude that vacuum is, indeed, the required 
base-geometry. 

As a final remark, it is also important to notice that 
both the impedance and the axion may be chosen to be 
spacetime dependent - a case that is not easily covered 
by TO, or geometrical optics. 



Rubilar [34[ (in the limit where sharp fronts propagate 
with no birefringence). One can confirm this statement 
by observing that the material in Sec. 8 of Ref. 34] is 
skewonless and obeys the closure relation (I105[) . Accord- 
ingly, the simple optical response (|106l) is retrieved, as 
explained in the manual [3]. One should finally notice 
that our findings could be relevant to the study of QED- 
induced spacetime structures [35j, |36[ . 
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V. CONCLUSIONS 

Starting from the classification of all skewonless, linear 
media by Schuller et al. Je}, we derived an intuitive rep- 
resentation for all non-birefringent materials: a metric 
medium core was employed, to which three "dispersion- 
relation preserving" corrections were added. 

This result provides a powerful inverse design tool. 
Once the optical metric of a material is known from prop- 
agation, the constitutive metric is fixed up to a constant: 
enforcing Q a @ cx G a @ is a very important step. Then, 
some freedom is still available in the corrections, namely 
an axion and two (anti-)selfdual bivectors (A and A). 

An explicit form for A and A was also developed, and 
was found to depend on the metric's signature. Remark- 
ably, the (anti-)selfdual bivectors were observed to van- 
ish in the Lorentzian (3,1) case. Accordingly, a reduction 
of all hyperbolic, skewonless, non-birefringent materials 
to metric media was achieved. The impact of this result 
both in prcmctric electrodynamics and in transformation 
optics was demonstrated ( Section II VI) . 

Considering further research, it would be interesting 
to study the medium tensor (|25[) more in detail. Indeed, 
the Fresnel equation (|2"5)l features interesting examples of 
separable (uniaxial) and non separable birefringent me- 
dia. In this context, one might even wonder whether the 
structure (|23|) covers all skewonless media. On a differ- 
ent note, the inclusion of a finite skewon in our analysis 
would also constitute a stimulating challenge. At the 
level of the Fresnel equation, a skewon term can be ac- 
counted for straightforwardly, thanks to the separation 
(D.2.42) in Ref. [3]. Additionally, a Segre classification 
of the full medium tensor still exists, even if x IJ cannot 
be diagonalized in general. Nonetheless, it is question- 
able that the ensuing classification remains manageable 
in a simple way. 

We believe that this work could spark some interest be- 
yond the linear-optics community. For instance, Eq. (|25p 
does encompass the nonlinear medium of Obukhov and 



Appendix A: Non-birefringent normal forms from 
Schuller's classification 



In this appendix, we classify all non-birefringent media 
by means typical representatives (normal matrix forms) . 
There exist in total five normal forms ("NBR" classes), 
which are obtained as limits from six specific "meta- 
classes" in Schuller's classification 0. As an aside, one 
should observe that some of the results reported here also 
appear in Section UlIDl 

The first NBR class is encoded by the matrix 





" — T 








a 





0" 







— T 








0" 













— T 








a 


Xi = 


a 








T 













a 








T 













a 








r 



(Al) 



which is obtained from metaclass I, with t\ = T2 = t% and 
o~i = o~2 — 03 . Equation (|A1[) covers all hyperbolic metric 
media, together with an axion term. The corresponding 
Fresnel equation is the free-space relation 



^2 



fc 3 



(A2) 



The second NBR class is attained from metaclass VII: 

A 6 



Xi; 



A 5 
A 3 





Ai 






A 6 

A 4 



0' 

A 4 

A 2 

A 5 
A 3 



A 2 





Ai. 



(A3) 



with the additional constraints 



Ai = si(A 2 - A 6 ) + s 2 A 5 , 
A3 = s 3 (A 4 - Ae) + S4A5 , 



(A4) 
(A5) 
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where si, s 2 , S3 and s 4 are signs. The associated disper- 
sion relation reads: 



S4 fc 3 2 ) 2 







(A6) 



so that the optical metric has signature (4,0) when S2 — 
S4 = +1, and (2,2) otherwise. From equation (| A3|) one 
can see that imposing A2 = A4 = \q reduces this NBR 
class to (4,0) or (2,2) metric media - with the addition 
of an axion piece. 

The remaining three NBR classes are "degenerate" , in 
that pairs of metaclasses give the same non-birefringent 
limit. This can be explained as follows: cr, and Tj are 
the real and imaginary parts of a complex eigenvalue Zi ; 
the symbol Ai denotes a real eigenvalue. Consider, then, 
some NBR matrix, which is derived by setting Tj = in 
a metaclass "A" . Clearly, enforcing Xj — in a metaclass 
"B" with one complex eigenvalue less, yields the same 
"degenerate" result. 

The third NBR class is given by 



(A7) 



This matrix is obtained from metaclass IV with (n — 
0, o-i = ±Ai + A 2 ) or (t 2 = 0, a 2 = ±Ai + A 2 ). If 
the second set of constraints is used, certain spacetime 
directions must be exchanged in order to retrieve Eq. 
(|A7[1 . Alternatively, the above normal form is the limit 
(Ai = 0, A 2 = ±A 3 + A 4 ) or (A 3 = 0, A 4 = ±A a + A 2 ) of 
metaclass VI. 

The fourth NBR class has the representation 



±Ai + A 2 













±Ai + 


A 2 


0" 





— r 








a 











Ai 








A 2 


±Ai + A 2 




















a 








r 











A 2 








Ai. 



Xiv 



A 3 

Ai 

±Ai + A 2 

A 4 

A 2 



A 4 

A 2 



A 3 

Ai 



(A8) 



which is achieved from metaclass VI with (n = 0, <T\ = 
±Ai + A 2 ), or from metaclass VII with (A5 = 0, X@ = 
±Ai + A 2 ). Other possibilities do exist; however, they 
correspond to a mere relabelling of coordinates. 
The fifth and last NBR class has the normal form 



(A9) 



where t\ = ±1. Eq. (IA9j) is attained from metaclass 
XVII by taking (r 2 = 0, A 2 = ±(A 3 — cr 2 )). Equivalently, 



"0 








Ai 




















A 5 











±(A 3 -A 5 ) 








A3 


A3 








ei 











As 




















Ai 








±(A 3 - A 5 ). 



one can start from metaclass XVIII, and set (A 4 = 0, 
A 2 = ±(A 3 - A 5 )) or (A 2 = 0, A 4 = ±(A 3 - A 5 )). 

The Fresnel equation of NBR classes iii, iv and v reads 



(±-fc 3 + fcifc 2 ) 2 = , 
c 



(A10) 

and corresponds to an optical metric with signature (2,2). 



Appendix B: The geometry of the NBR classes 

We shall now describe how to derive the general for- 
mula (|25[) from the NBR matrices of Appendix \K\ Given 

the dispersion relation of a normal form, the metric Q is 
fixed up to a constant. This specifies the "core" of the 
medium response: a "remainder" (or "residue") is left 
out, which is simplified - fully, or to some extent - by a 
suitable choice of axion. At this point, the problem for 



NBR class i solved immediately, as one finds A 
For NBR class ii, the residue reads: 



= A = 0. 



"0 














■ 





s a2 a 2 








o 2 



















ai 























ft 2 








s a2 a 2 











Ql 








s ai ai_ 



(Bl) 



where s Q2 and s ai are signs and satisfy s a2 = s 2 s 4 s Ql . 
By inspection, one can conclude that the structure (|B1[) 
does indeed consist of two bivectors A and A. 

For NBR classes iii, iv and v, the remainder is of the 
form (in the case of NBR class v, one must re-label the 
coordinates) : 



Ax, 



i/iv/v 



0" 

a 7 





7 f3 





(B2) 



Again, a simple calculation reveals that this matrix can 
be represented by A and A (cf. Sec. lIIIDj) . Consequently, 
Eq. (|25p must encompass all non-birefringent media. 

As a final remark, explicit expressions for the bivectors 
can be deduced from the residues in this Appendix. The 
resulting quantities (A, A and Q) always comply with 
the conditions (|2^)) -(|5T |) . This constitutes an excellent 
consistency check. 



Appendix C: Schuller's classification and premetric 
electrodynamics 

The classification proposed by Schuller et al. Q is 
(background) metric free by construction. Yet, it relies 
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on a linear medium response Q^ va " which is a true ten- 
sor: 

Qfi'u'a'ff' = l^l^l^l^ QM"«/3 . ( C1 ) 

In the premetric formalism of Sec. IIII1 one must derive 
Qvap f rom ^he can onical x™ ■ For this purpose, Ref. 
§ develops a scheme to eliminate the "extra" density fac- 
tor in Eq. (|103l) . The present appendix analyzes this pro- 
cedure; it demonstrates its consistency, but reveals that a 
specific uniqueness proof is still missing. Equivalently, a 
(cautiously) positive answer is found to the question: can 
Schuller's classification be employed in premetric electro- 
dynamics? 

Schullcr ct al. examine the quantity | det(r2 7J )| 1 / 6 , and 
show that it is a scalar density: 

| det(n I ' J ')\ 1 ^ = I det(L^)| +1 | det(f* /J )| 1/6 , (C2) 

where det(fi 7J ) is a 6x6 determinant. As an aside, it 
is interesting to make a comparison with the true scalar 
I det(XpM)|, satisfying: 

|det(x£f)l =\det(xl J u )\ . (C3) 

Resuming the derivation of Ref. Q , the weight of — 1 ob- 
tained in (IC2|) is used to compensate the w = +1 scaling 
of e a ^ 5 (see Eq. (fTUTT) ): 

v ap 7 s ._ | det^- 7 )! 1 ^/?^ ( C4 ) 



Remarkably, a Levi-Civita-type tensor, ii"^ 7 *, is there- 
fore achieved, which, under a change of frame, obeys 

A similar procedure is exploited to link to the 

canonical Xpm"^ % so as to attain the relation 



|det( X ^)|V6 |det(fi")|V6 ' v ; 

Here, the expression on the right hand side can be found 
in Q , while the expression on the left hand side is a slight 
generalization comprising those Xpm*' 3 that are not uni- 
modular in a 6x6 sense. Thus, combining the transfor- 
mation properties (|103IC1IC2IC3I) . one reaches the im- 
portant conclusion that the density-elimination scheme 
proposed by Schuller et al. is fully consistent. 

Besides this positive outcome, one can observe that 
the two fractions in (|C6[) are functions of two separate 
objects, tl^ap a nd x™"^ ■ Consequently, one can never 
uniquely derive a tensorial fl^"^ from the more fun- 
damental Xpm"' 3 . Accordingly, a density- free formula- 
tion, even if consistently constructed, does not appear 
to emerge unambiguously from the axiomatic structure 
of pre-metric electrodynamics Q - at least so far. 

In this paper, there exists an intuitive way to account 
for densities: once a metric Q a P is extracted from f2^ l/Q!/3 , 
the identity 

^\det(Q-^,)\ = |det(L";)|- 1 v /|det(Q^)| , (C7) 

can be used. Consequently, one can readily make a tran- 
sition to Xpm"^, as demonstrated in Eq. (I104|) : the prob- 
lem of uniqueness found in (IC6[) is accounted by the pref- 
actor M, and does not affect our results fSec. IIV A")) . 
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